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We analyze a novel method that uses fixed, minimal physical resources to achieve generation 
and nested purification of quantum entanglement for quantum communication over arbitrarily long 
distances, and discuss its implementation using realistic photon emitters and photonic channels. 
In this method, we use single photon emitters with two internal degrees of freedom formed by an 
electron spin and a nuclear spin to build intermediate nodes in a quantum channel. State-selective 
fluorescence is used for probabilistic entanglement generation between electron spins in adjacent 
nodes. We analyze in detail several approaches which are applicable to realistic, homogeneously 
broadened single photon emitters. Furthermore, the coupled electron and nuclear spins can be used 
to efficiently implement entanglement swapping and purification. We show that these techniques 
can be combined to generate high-fidelity entanglement over arbitrarily long distances. We present a 
specific protocol that functions in polynomial time and tolerates percent-level errors in entanglement 
fidelity and local operations. The scheme has the lowest requirements on physical resources of any 
current scheme for fully fault-tolerant quantum repeaters. 

PACS numbers: 03.67.Hk, 03.67.Mn, 78.67.Hc 



I. INTRODUCTION 

Quantum communication holds promise for transmit- 
ting secure messages via quantum cryptography, and for 
distributing quantum information However, expo- 
nential attenuation in optical fibers fundamentally lim- 
its the range of direct quantum communication tech- 
niques , and extending them to long distances remains 
a conceptual and technological challenge. 

In principle, the limit set by photon losses can be over- 
come by introducing intermediate quantum nodes and 
utilizing a so-called quantum repeater protocol Such 
a repeater creates quantum entanglement over long dis- 
tances by building a backbone of entangled pairs between 
closely-spaced quantum nodes. Performing an entangle- 
ment swap at each intermediate node 4] leaves the outer 
two nodes entangled, and this long-distance entangle- 
ment can be used to teleport quantum information [^|^ 
or transmit secret messages via quantum key distribu- 
tion 0. Even though quantum operations are subject to 
errors, by incorporating entanglement purification 0,0] 
at each step, one can extend entanglement generation to 
arbitrary distances without loss of fidelity in a time that 
scales polynomially with distance This should be 

compared to direct communication, which scales expo- 
nentially, making it impractical for long distances. 

Several approaches for physical implementation of a 
quantum repeater protocol have been proposed. Early 
work was based on systems of several atoms trapped in 
high finesse optical cavities 0, Q • Such systems 
can form a quantum network with several quantum bits 
(qubits) per node, and are particularly suitable for effi- 
cient implementation of the pioneering proposal of Ref. 
U . In this approach, quantum communication over thou- 



sand kilometer distances requires seven quantum bits per 
node, which must be coherently coupled to perform local 
quantum logic operations, i.e. a seven qubit quantum 
computer. The specific implementation of these early 
ideas involved the techniques of cavity QED for interfac- 
ing stationary and photonic qubits and for performing 
the necessary quantum logic operations [l3l . Recent 
related work pointed out that long-distance entanglement 
can be implemented via probabilistic techniques without 
the use of ultra- high finesse cavities |^ 0| , while local 
operations can be accomplished via short-range interac- 
tions involving e.g. interacting trapped ions. However, 
few-qubit registers are still technically very difficult to 
construct, and the difficulty increases drastically with the 
number of qubits involved. At the same time, a novel ap- 
proach based on photon storage in atomic ensembles [15| 
and probabilistic entanglement is also being actively ex- 
plored. In comparison with systems based on many-qubit 
nodes, this approach offers less error tolerance and re- 
quires a longer communication time. Realization of a 
robust, practical system that can tolerate all expected 
errors remains therefore a challenging task. 

In a recent paper 1^ we proposed a quantum repeater 
protocol which could be implemented using the electronic 
and nuclear degrees of freedom in single photon emitters. 
Here we present further details of the proposal described 
m Ref. [ig, and we compare our methods to alternative 
strategies. We show that our repeater protocol requires 
only two effective quantum bits at each node. This is the 
minimum requirement on physical resources which still 
allows active error correction. As a specific implemen- 
tation, we consider nodes formed by a single quantum 
emitter with two internal degrees of freedom. A pair of 
electronic spin sublevels allows for state-selective optical 
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excitation (see inset in Figure la), and a proximal nuclear 
spin provides an auxiliary memory. State-selective fluo- 
rescence is used for probabilistic entanglement generation 
between electronic spin degrees of freedom. We analyze 
in detail and compare several approaches for probabilis- 
tic entanglement generation, focussing on the feasibility 
of their implementation using realistic photon emitters. 
Once electronic spin entanglement is generated, the cou- 
pled electron and nuclear spin at each node can be used 
to efficiently implement entanglement swapping and pu- 
rification. We show that these techniques can be com- 
bined to generate high-fidelity entanglement over arbi- 
trarily long distances. We present a specific protocol that 
functions in polynomial time and tolerates percent-level 
errors in entanglement fidelity and local operations. 

Our approach is stimulated by recent experimental 
progress in single photon generation by means of single 
quantum emitters, including atoms and ions as well as 
impurities and nanostructures in solid state devices. Al- 
though our approach is relevant to atomic systems, such 
as single atoms trapped in a cavity |17| or single trapped 
ions |l2l |. it is particularly suitable for implementation 
with solid-state emitters, for ex amp le impurity color cen- 
ters 0,01 and quantum dots |2nLl2ll |. These devices of- 
fer many attractive features including optically accessible 
electronic and nuclear spin degrees of freedom, potential 
opto-electronic integrability, and fast operation. 

The paper is organized as follows. First, we will discuss 
techniques for entanglement generation. For clarity, we 
will present our results in the context of nitrogen- vacancy 
(NV) centers in diamond, and discuss alternative im- 
plementations at the end. Realistic imperfections, such 
as homogeneous broadening and limited selection rules, 
motivate a novel entanglement generation scheme based 
on state-selective Rayleigh scattering and interferometry. 
We calculate the success probability and entanglement 
fidelity for this scheme as implemented in NV centers, 
and compare this scheme to alternative schemes based on 
Raman scattering or optical tt pulses, with success condi- 
tioned on detection of one or two photons. Next, we will 
show how hyperfine coupling between the electron spin 
and proximal nuclei permits entanglement swapping and 
purification. Performing these operations in parallel and 
employing a nesting scheme, we calculate the fidelity ob- 
tained and the time required to generate it as a function 
of distance. In addition, we compare this scheme to pi- 
oneering proposals [3, 0| for fault-tolerant quantum re- 
peaters. Finally, we quantitatively discuss the feasibility 
of implementing a quantum repeater using NV centers, 
and elucidate alternative material systems which satisfy 
the requirements of our technique. 



II. ENTANGLEMENT GENERATION 

The initial step in our scheme is entanglement gener- 
ation between the electron spins of two emitters sepa- 
rated by a distance Lq. In principle, entanglement can 



be generated probabilistically by a variety of means, e.g., 
Raman scattering [l^ |23, '^^ polarization-dependent 
fiuorescence 0. However, for our repeater protocol it 
is essential that the optical transition be independent of 
the nuclear spin state, and solid state emitters do not al- 
ways allow Raman scattering or polarization-dependent 
fiuorescence which fulfills this requirement. We there- 
fore consider an entanglement mechanism based on state- 
selective elastic light scattering as shown in Figure 1. 
Elastic light scattering places few restrictions on selection 
rules, and permits nuclear-spin-independent fiuorescence 
as we discuss below. 



A. Properties of single color centers 

Our entanglement generation scheme is applicable to 
a wide variety of physical systems, requiring only the 
simple level structure illustrated in Fig. For clarity, 
we will present it first using a concrete example: the 
nitrogen-vacancy (NV) center in diamond, which has the 
specific level structure shown in Fig. [21 This example 
illustrates many generic features common to other solid- 
state emitters. 

NV centers represent a promising physical system be- 
cause of their strong optical transition around 637 nm 
and optically accessible electron spin states. In partic- 
ular, the ground state {Ai symmetry class of the C^y 
group) has three electronic spin-states which we label 
I — 1), |0) and |1) according to their spin angular mo- 
mentum along the symmetry axis of the crystal (Mg). 
Spin-orbit and spin-spin effects lead to a splitting of |0) 
from I ± 1) by 2.88 GHz. Since we only require two elec- 
tronic spin states, |0) and we isolate these two states 
from the | — 1) state by either adding a small magnetic 
field to shift the energy of the | ± 1) state, or by us- 
ing appropriately polarised ESR-pulses. As spin-orbit 
and spin-spin effects are substantially different for the 
optically excited state (E symmetry class), the strong 
transition from the Mg = sublevel of the ground or- 
bital state can be excited independently of the other Mg 
states. Although there is evidence for photo-bleaching at 
low temperatures, current models indicate that crossover 
into the dark metastable state occurs primarily from the 
Ms ~ ±1 excited states |23|- Furthermore, crossover 
into the trapping state is a detectable error. In the re- 
peater protocol described below we perform a series of 
measurements on the electronic spin. During these mea- 
surements, the dark state will not produce any fluores- 
cence, revealing the error. Shelving into the metastable 
state will thus influence the time (see Appendix) but not 
the fldelity associated with the repeater protocol. Con- 
sequently, we assume that we are only near resonance 
with a single state |e) which has Ms = 0, and neglect 
photo-bleaching effects. 

The electron spin degree of freedom suffices to generate 
entanglement between adjacent NV centers. To propa- 
gate entanglement to longer distances, we will make use 
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FIG. 1: (a) Generic level structure showing the state-selective 
optical transitions and electronic spin sublevels required for 
entanglement generation, (b) Setup used to create entangle- 
ment. The two emitters act as state dependent mirrors in an 
interferometer. The outputs of the cavities (a\ and 02) are 
combined on a beamsplitter. By proper alignment of the in- 
terferometer the photons always exit through the (ai + a2)/\/2 
port if both centers are in the scattering state |0). A detection 
of a photon in the (ai — a2)/\/2 mode thus leads to an entan- 
gled state, (c) Scheme for balancing the interferometer. Each 
node is optically excited by a laser pulse which first reflects 
off the other node, so that the optical path lengths for the two 
excitation/ emission paths are identical. 



of an auxiliary nuclear degree of freedom {||) , ||)} which 
will be used for storage of quantum information during 
the repeater protocol. In NV centers, this nuclear degree 
of freedom can arise from a nearby carbon-13 impurity or 
directly from the nitrogen-14 atom that forms the cen- 
ter. The large energy separation between the |0) and 
I ± 1) states exceeds the hyperfine interaction by an or- 
der of magnitude, decoupling the nuclear and electronic 
spins. The energy levels can thus be described by prod- 
uct states of the two degrees of freedom. Furthermore, 
in states with Mg = 0, the energy is independent of the 
nuclear state. Finally, a small magnetic field ^ 10 — 100 
Gauss allows spectral resolution of the Mg = ±1 states 
without producing significant nuclear Zeeman splitting. 
The optical transition between |0) and |e) is thus disen- 



tangled from the nuclear spin state. Consequently, the 
nuclear spin can be used to store entanglement while the 
|0) — |e) transition is used to generate another entangled 
pair of electron spins. 
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FIG. 2: The relevant electronic and nuclear states of the cou- 
pled NV center and C impurity nuclear spin. The electron 
spin states can be coupled by ESR microwave fields near 2.88 
GHz, while the nuclear spin states can be addressed by NMR 
pulses on the 130 MHz hyperfine transition. A laser applied on 
resonance with the Ma = optical transition produces strong 
fluorescence for the entanglement scheme; it will also weakly 
excite the Ms = ±1 transitions. In our model we assume that 
only the Ms = ±1 states decay to the shelving state \W) at 
rate Ts- 



B. Entanglement protocol 

To implement the entanglement scheme, each NV cen- 
ter is placed inside a photonic cavity, whose output is 
coupled to a single-mode photonic fiber (note, however, 
that cavities are not essential for this proposal, see be- 
low). Fibers from adjacent NV centers enter two ports 
of a beamsplitter, and entangled electron spin pairs are 
obtained conditional on detection of a photon in an out- 
going arm of the beamsplitter. 

Specifically, our protocol for entanglement generation 
relies on scattering on the closed optical transition be- 
tween |0) and |e). This scattering does not change the 
state of the NV center; the centers essentially act as mir- 
rors reflecting the light only if they are in the state |0). 
We assume that each of the centers is initially prepared 
in the state (|0) + |l))/\/2, so that the total state is 
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Since there is no light scattering from state we can 
exclude the component if we detect any scattered 
photons. In state |00), both centers act as mirrors, so 
that by balancing the interferometer in Fig. ^ (b) we can 
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arrange for the photons to leave through a specific port 
. A photon detection in the opposite port D_ can 
thus only arise from the |01) and |10) states and produces 
an entangled superposition of these two states. 

Balancing and stabilizing an interferometer over tens 
of kilometers as required for the implementation of this 
protocol represents a considerable challenge. Using a 
method analogous to the plug- and-play system used in 
quantum key distribution [23 , we can reduce this require- 
ment to stabilization of a small interferometer locally at 
each detector. Suppose that we wish to generate entan- 
glement between repeater nodes Ri and i?2- Employing 
fast optical switches, we excite i?i by sending a pulse 
of light toward i?2, where the light is reflected and sent 
past the detector to Light emitted from i?i follows 
the same path back to the detector. Similarly, we excite 
i?2 by sending a pulse of light toward a fast switch at Ri . 
The two paths thus cover the same path length between 
the nodes and we are insensitive to fluctuations in the 
path lengths as long as these fluctuations happen on a 
time scale which is longer than the time it takes for the 
pulses to travel between the stations. Alternatively one 
could change to a protocol which relies on the detection 
of two photons instead of one. In such protocols the sen- 
sitivity to changes in the path lengths can be reduced 
considerably [2al27Ll2^. 



C. Entanglement fidelity in the presence of 
homogeneous broadening 

We now describe this process mathematically, calcu- 
lating the fidelity of the entangled pair produced by our 
protocol, as well as the probability for it to succeed. Our 
analysis incorporates dominant sources of error in solid- 
state systems; in particular, we account for effects of ho- 
mogeneous broadening on the optical transition. 

Our model assumes that the NV centers are excited by 
a weak driving field applied between the states |0) and |e) 
with Rabi frequency 17; the excited states then radiatively 
decays at a rate 7. To describe the effect of homogeneous 
broadening on the optical transition we assume that the 
energy of the excited level fluctuates with a characteris- 
tic time which is slow compared to the optical frequency 
and much shorter than the lifetime of the excited state. 
In this approximation the broadening can be described 
by including a time-dependent detuning A(i) with white- 
noise characteristics: (A(t)) = 0, (A(i)A(t')> ^T5{t-t'). 
Below we shall be working in the limit of weak driving 

^ 7 -f- r. In this limit the light emitted from a center 
consists of two contributions: (i) a coherent part cen- 
tered around the frequency of the driving laser with a 
width given by the width of the driving laser, and (ii) 
an incoherent part centered around the frequency of the 
transition with a frequency width of 7 -I- F. The relative 
weight of these two contributions is 7 : F. With consid- 
erable broadening of the optical transition F > 7 it is 
therefore essential to filter out the incoherent scattered 



light with a frequency filter to get a high fidelity. To 
filter out the incoherent light and obtain a high collec- 
tion efficiency we assume that the centers interact with 
an optical cavity with a coupling constant g and a de- 
cay rate k. We emphasize, however, that good cavities 
are not essential for our proposal: we only require suffi- 
cient collection efficiency and frequency selectivity, which 
could also be obtained by collecting the light with a lens 
and sending it through a narrow frequency filter. In gen- 
eral the weak drive may be detuned from the excited 
state, which would simplify the filtering of coherent from 
incoherent light. However, off-resonant excitation would 
require a stronger driving field, making it harder to avoid 
stray light reaching the detectors. For simplicity we only 
discuss the situation where, on average, the driving field 
and cavity mode are resonant with the center. 

The combined NV-center cavity system is then de- 
scribed by the Hamiltonian 

H ^ A(t)|e)(e| + ^(|0)(e| -f |e)(0|) + gct|0)(e| +.g|e)(0|c, 

(2) 

where c is the photon annihilation operator for the field in 
the cavity. In the Heisenberg picture, decay terms can be 
included by considering the quantum Langevin equations 
of motion for the atomic operators (Ty — \i) 



dc 
dt 



dt = [-i-^m-o. 



M + — ) ((Tee - O-QO) + ^Oe 



(3) 



(4) 



= -l^ee + + gC^)^Oe + fi-C. ) + Fee (5) 



where the noise Fc is the incoming vacuum noise leading 
to cavity decay at rate k and the other noise operators 
Foe , Fee represent the effect of other optical modes that 
lead to decay. 

We obtain an appropriate solution to the quantum 
Langevin equations by noting that, in the limit of weak 
driving, f2 7, there is virtually no population of 
the excited state, ctoo — <Tee ~ (Too- The solution can 
then be written in the form c = aaoo+noise and aoe = 
/?(Too+noise, but the equations for a and P are compli- 
cated due to the noise A{t). By averaging the Langevin 
equations over the noise one can, however, find simple 
equations for various moments of a and (3, and by taking 
steady state solutions of the averaged equations we find 



-2gn 



K(7 + F)(l + 4g2/'^(7 + r)) 



K2(7+r)2 



^+K(7+r)yv^ (7+r)(7+K) + K(7+r) 

-in 



(7 + F)(l + 4gV.j(7 + F))- 



(6) 



(7) 



(8) 
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Note that in the presence of ho moge neous broadening 
r ^ 0, the moments do not factor, |ap 7^ signifying 
incoherent scattering of hght into the cavity. 

We now apply the entanglement generation protocol, 
and use our mathematical model to predict the aver- 
age density matrix components of the NV center elec- 
tron spins. In our scheme, we combine the output of 
the two cavities on a beamsplitter and select the de- 
sired entangled state by conditioning on a click in detec- 
tor £)_, described by the photon annihilation operator 
d- — y%K/2(ci — £2). Here, subscripts one and two re- 
fer to the two NV-centers we are trying to entangle, ( is 
the total collection and detection efficiency for photons 
leaving the cavity, and we have omitted the contribution 
from vacuum noise operators. To describe the effect of 
the detection, we use the quantum jump formalism(29l|. If 
the system starts out in state |\I'imt), the density matrix 
element pij at time t can be found by 

P^At) = {■^nut\d-{ty\j){t\d^{minit)St/SP, (9) 

where the time argument t is included to emphasise the 
time dependent Heisenberg operators, and where SP is 
the probability to have a click during a time dt, 



SP = {^init\d-{t)^dAminit)St. 



(10) 



Our entanglement generation scheme relies on interfer- 
ence to eliminate D_ detection events coming from the 
initial state |00). However, according to our formalism, 
if we start out in an initial state |00) the probability to 
have a click is given by 



SP (5t(00|dLd_|00) 



KCSti |a|2 - 



(11) 



where we assume the noise is independent for the two 
centers. This expression vanishes only for coherent scat- 
tering of Hght into the cavity, i.e. |aP = |ap or F = 0. In 
the presence of broadening, there is a finite probability 
that light will be detected from the |00) state. Similarly, 
r > leads to a finite probability for incoherent scat- 
tering from |01) and |10). Homogeneous broadening thus 
reduces the fideUty {F = (*idoai|/3|*idoai), (where l^'idoai) 
denotes the ideal entangled state) by 



F = 



lap 



1 



2 7 + r 7 + K 1 452/^(7 -I- r) ' 



(12) 



Here we are interested in the limit where the fidelity is 
close to unity and we shall therefore assume |aP ~ jap 
in the calculation of other noise sources below. 



D. Other errors 

In addition to the error caused by homogeneous broad- 
ening, there is also a reduction in fidelity caused by mul- 



tiple spontaneous emission events from the centers. This 
fidelity can conveniently be expressed in terms of the to- 
tal emission probability 

(7 + r)[i + v/«:(7 + r)]' ^'""^ 

where to is the duration of the applied laser pulse. In the 
absence of homogeneous broadening, multiple excitations 
result in a fidelity 

1 c-^™ (1-^/2) , , 

^ (14) 



^^2 
and success probability 



P - (1 - exp(-ePem/2))/2. 



(15) 



The total collection efficiency can be expressed as e = 
C^cav with the probability to emit into the cavity given 

by 



P — 



(16) 



n- 452/^(7 + r)' 

This treatment has neglected the possibility of distin- 
guishing multiple photon detection events. If our detec- 
tor can resolve photon number, we can use the informa- 
tion to improve our protocol. In particular, a detection 
in the mode described by d+ = y(>«/2(ci -I- C2)-f noise 
has no effect on the component of Eq. that we are 
interested in, since d+(|01) + |10)) oc (|01) + |10)). Fur- 
thermore, a detection in this plus mode contains contri- 
butions from 1 00), so it yields no useful information. On 
the other hand, detection events in the mode described 
by d- change the sign of the superposition state, since 
d_(|01) + |10)) cx (|01) - |10)), and d_(|01) - |10)) cx 
(|01) -I- |10)). Consequently, the optimal strategy is to 
change the phase of the entangled state when an even 
number of photons is detected. The resulting fidelity is 

-Pe™(l-e) 

• (17) 
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Finally, we must include the effect of two other sources 
of noise: dark counts and electron spin dephasing. In 
the limit of small success probability P the dark count 
introduces an incoherent admixture of the initial state 
and thus leads to a reduction in fidelity Pdark/P where 
-Pdark = 7dc^o is the dark count probability. Electron spin 
dephasing makes the state decay towards a state with fi- 
delity 1/2 at a rate 27e, yielding a reduction in the fidelity 
of 7e times the total time of the experiments. Typically, 
this total time will be dominated by the classical com- 
munication time between nodes, tc- 

Putting these considerations together, we find that 
the entanglement scheme succeeds with probability P = 



(1/2) (1 



e-fem/4, producing the state |'I'_) 



in time Tq « (to + tc)/ P with fidelity 

Fo = i(l + e-^™(i-)) -7e(to + <c) 



to 



3 F 



1 



-Idc p 2F + 7K + 7I+. 



V«(7 + r) 



(18) 
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For realistic emitters placed into a cavity with either 
a narrow linewidth k <C 7 or a large Purcell factor 
4(7^/(k(7 + r)) 3> 1, the first two terms should domi- 
nate the error. 



III. COMPARISON TO OTHER 
ENTANGLEMENT GENERATION SCHEMES 

The entanglement generation scheme that we have pre- 
sented so far is the scheme that we believe to be best 
suited to NV centers. For other systems, this may not 
be the case. In particular, the presented scheme has 
two primary drawbacks: (1) it relies on resonant scatter- 
ing, making it difficult to filter fluorescence photons from 
the applied laser field; (2) to avoid loss of fidelity from 
incoherent scattering, one must detect only a narrow 
frequency interval in the scattered light. Other entan- 
glement methods present different problems which may 
prove easier to resolve or other methods may be better 
suited for different physical systems. Consequently, we 
now briefly compare the resonant scattering scheme pre- 
sented above to alternate techniques. 

A. Raman transitions 

One of the first schemes considered for probabilistic 
entanglement generation [isl 123 . |2^ used Raman transi- 
tions in three level atoms. In such schemes, an electron 
spin flip between non-degenerate ground states |0) and 
|1) is associated with absorption of a laser photon and 
emission of a frequency shifted Raman photon. After 
interfering the emission from two atoms, detection of a 
Raman photon projects the two-atom state onto a state 
sharing at least one flipped spin. To avoid the possibility 
that both atoms emitted a Raman photon, the emission 
probability must be quite small Pem <C 1. In this limit, 
a photon detection event in detector D± results in an 
entangled spin state |^±). 

The Raman scheme can be implemented using either 
a weak drive between states |1) and |e) or with a short 
strong pulse which puts a small fraction of the population 
into |e). Since the latter is equivalent to the single de- 
tection TT-pulse scheme discussed below, here we consider 
only weak driving. The system can now be treated using 
the quantum Langevin-quantum jump approach formu- 
lated above. As before, homogeneous broadening on the 
optical transition leads to an incoherent contribution to 
the Raman scattered light, which reduces the entangle- 
ment fidelity in a manner similar to Eq. H12|) . Again, 
for optimal fidelity the coherent part should be isolated 
with a narrow frequency filter. If we assume perfect fil- 
tering and a small collection efficiency e ^ 1 the fidelity 
conditioned on a click is given by 

P = 1 - Pcm (19) 

with success probability P — PcmC- 



In the limit of large fidelity P « 1 , the Raman scheme 
has a success probability which is a factor of 4 higher 
than for our interferometric scheme. Furthermore, the 
Raman scheme has the advantage that stray light may 
be spectrally filtered from the Raman photons. Never- 
theless, because of the hyperfine interaction in state 
the transition frequency from |1) to |0) depends on the 
nuclear spin state. The associated detrimental effect on 
the nuclear coherence could potentially be avoided by us- 
ing simultaneous transitions from |1)| f) and | — 1)| i), 
which are degenerate. To our knowledge, however, fluo- 
rescence between |e) and |1) has not be observed, making 
it uncertain whether the Raman scheme can be imple- 
mented for the NV-centers [49| . 



B. TT pulses 

Time-gated detection offers an alternate method for 
distinguishing scattered photons from stray incident 
light. If an atom or NV center is excited by a suffi- 
ciently short, strong laser pulse, its population is coher- 
ently driven into the excited state. The excited state |e) 
then decays on a time scale I/7. When the decay time 
is much longer than the incident pulse length, the exci- 
tation light and the photon emitted from the atom are 
separated in time, and can thus be distinguished. Entan- 
glement is then generated conditional on the detection of 
one or two photons, as elucidated below. 



1. Single detection 

One particularly simple method for generating entan- 
glement using TT-pulses begins with each atom in a state 

cos(</))|l) +sin((/))|0). (20) 

The incident 7r-pulse excites the optically active state |0) 
to |e) with unit probability, and the spontaneously emit- 
ted photons are interfered on a beamsplitter and sub- 
sequently measured (as for the Raman scheme above). 
Provided (/) <C 1 we can ignore the possibility that both 
atoms are in state |0) . A photon detection in D± excludes 
the state preparing the system in |^'±). 

As with other entanglement schemes, high-fidelity en- 
tanglement generation requires filtering the incoherent 
scattering caused by homogeneous broadening of the op- 
tical transition. In previous sections, we have proposed 
to use a frequency fllter to separate the narrow peak (in 
frequency) of coherent scattered light from the broad in- 
coherent background. In the present case, filtering can 
be done in the time domain. In the excitation process a 
coherence is established between |1) and |e), and follow- 
ing the excitation this coherence (the off-diagonal density 
matrix element) decays at a rate r/2. By only condi- 
tioning on photons emitted a very short time after the 
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excitation, during which the coherence has not had time 
to decay, a high quahty entangled pair is produced. 

To describe this process mathematically, we again as- 
sume that the atom is placed inside an optical cavity. 
In contrast to our previous calculations, we assume that 
the cavity has a broad linewidth to ensure that gener- 
ated photons leave the system as fast as possible. In the 
limit g, "f, and T we can adiabatically eliminate the 
cavity by setting dc/dt — in Eq. Q so that if we omit 
the noise we obtain 



c{t) = aoe 

K 



(21) 



Inserting this expression into Eqs. I@J and we find 

&^e{t) - &oe{t = 0) exp (- ^"''^^ t] (22) 



(23) 



O-ee(t) = aee(t = 0)exp(-7cfft), 



where the effective decay rate jeS is the decay rate en- 
hanced by the Purcell effect 



7cff 



492 

7 1 + — 

K7 



(24) 



To find the fidelity of the entangled state created with 
this method we again use Eqs. © and (|10|l . For simplic- 
ity we only work in the limit of small collection efficiency 
e <C 1. Conditioned on a click at time t after the excita- 
tion, the fidelity of the entangled state is 



^^^cos2(0)(i + V" 



(25) 



and the probability to have a click during the short time 
interval from t to t + St is 



SP = 2e7eff<5isin2((/))e' 



-7cttt 



(26) 



where the collection efficiency e = C^cav is again given by 
the collection efficiency for the light leaving the cavity 
and the probability to emit into the cavity is now given 
by 



4g2/K7 



I 



Vk7 



(27) 



The success probability for a given fidelity now de- 
pends on the ratio between the broadening and the ef- 
fective decay rate r/7off. For F = 0, the procedure 
of initially transferring population from |1) to |0) and 
then applying a 7r-pulse between |0) and |e) is equivalent 
to a Raman transition, and Eqs. (|25I26I) indeed repro- 
duces the same relation between success probability and 
fidelity given in Eq. H19|l . In the limit of small broaden- 
ing, F <C 7cff, the TT-pulse scheme is advantageous over 
the interferometric scheme presented first. In particu- 
lar, for a fixed fidelity F w 1 the success probability is a 
factor of 4 higher. 



In the presence of broadening, however, the situation 
is different. To obtain a high fidelity we should detect 
only photons emitted within a short time T following the 
excitation. The average fidelity will then depend on two 
parameters and T. By optimizing these two parameters 
we find that for _F w 1 the fidelity is 



F = 1- 



F 



87cff 



(28) 



Since previous expressions H14|) and (|I9|I for 1 — F de- 
pended linearly on P, this represents a much faster de- 
crease in the fidelity. The tt— pulse scheme is thus less 
attractive for homogeneously broadened emitters. 



2. Double detection 

If the collection efficiency is very high, it may be an 
advantage to rely on the detection of two photons instead 
of one 26, 27„,3Qei31J. In this scheme, both atoms are ini- 
tially prepared in (|0) -I- |l))/-\/2 and a 7r-pulse is appHed 
between |0) and |e). Following a detection in D± the pop- 
ulations in states |0) and |1) are interchanged and another 
TT-pulse is applied between |0) and |e). Conditioned on 
clicks following both tt pulses we can exclude the possi- 
bility that the atoms were initially in the same state and 
we are left with |5'±) conditioned on appropriate detector 
clicks. In the absence of homogeneous broadening, this 
protocol produces an entangled state with fidelity F=I 
with probability P = 6^/2. The double-detection scheme 
thus avoids the multiple photon emission errors inherent 
in the single-detection schemes. 

With broadening of the optical transition, this is no 
longer the case. For F « f , the relation between fidelity 
and success probability is now given by 



F = 1 - 



2P. 



7cffe 



(29) 



Again the fidelity decreases more rapidly with the success 
probability than for the Raman and resonant scattering 
scheme, making it less useful for our purpose. 



C. Summary 

The best choice of scheme depends on the specific phys- 
ical situation. The two 7r-pulse schemes are advantageous 
if the broadening is negligible. In particular, in the limit 
where we can ignore all errors except the photon attenu- 
ation, the double detection scheme results in the highest 
fidelity entangled pair. With low collection efficiency or 
large distances between emitters, the double detection 
will have a very small success probability because of the 

factor, and it may be advantageous to rely on a single 
detection scheme. 
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The TT-pulse schemes are less attractive if we are hmited 
by homogeneous broadening of the optical transition be- 
cause the fidelity decreases rapidly with the success prob- 
ability. Better results are obtained for the resonant scat- 
tering or Raman schemes. When possible, the Raman 
scheme offers the best solution. The frequency-shifted 
Raman scattering allows frequency filtering of the incom- 
ing light; in addition the success probability is four times 
higher than for the resonant scattering scheme. But, as 
mentioned above, it is not always possible to drive Raman 
transitions , and it may be hard to achieve Raman tran- 
sition which are independent of the nuclear spin state. 
For this reason we believe that the resonant scattering 
scheme is most promising in the particular case of NV- 
centers. 

Finally we wish to add that the calculations we have 
performed here assume a specific model for the broad- 
ening (short correlation time for the noise). With other 
broadening mechanisms, e.g. slowly varying noise, these 
considerations will be different. 



L„ L„ < L„ ► 



>• 




t T 

Measure in Bell state basis 
of each two-qubit system 

FIG. 3: Entanglement propagation by swapping. To gener- 
ate an entangled nuclear spin pair (black) over the distance 
Ln+i ~ 2Ln + Lo, we first generate nuclear entanglement 
over the first and second pair of repeater stations which are 
distance L„ apart. The electron spin (red) is then used to 
generate entanglement between the middle stations, separated 
by distance Lq. Nuclear (electron) spin entanglement is illus- 
trated by solid (dashed) lines. Entanglement swapping is per- 
formed by measuring in the Bell state basis of the two-qubit 
system. 



IV. ENTANGLEMENT SWAPPING AND 
PURIFICATION 

Using one of the procedures outlined above, electron 
spin entanglement can be generated between adjacent 
pairs of nodes. We now discuss a means to extend the 
entanglement to longer distances. 

A. Swapping 

After entangling nearest-neighbor electron spins, the 
electron spin state is mapped onto the auxiliary nu- 
clear spin qubit for long-term storage using the hyper- 
fine interaction. This operation leaves the electronic 
degree of freedom available to generate entanglement 
between unconnected nodes, as illustrated in Figure |31 
By combining optical detection of individual electron 
spin states jS^I and effective two-qubit operations as- 
sociated with hyperfine coupling of electronic and nu- 
clear spins Js^, we may projectively measure all four Bell 
states in the electronic/nuclear manifold associated with 
each emitter. The outcomes of the Bell state measure- 
ments reveal the appropriate local rotations to obtain a 
singlet state in the remaining pair of nuclear spins, im- 
plementing a deterministic entanglement swap 0, . By 
performing this procedure in parallel, and iterating the 
process for N cx log2{L/Lo) layers, we obtain the de- 
sired nuclear spin entanglement over distance i in a time 
cx Llog2 {L/Lo). 



of the repeater scheme. By performing local operations 
and measurements, it is possible to distill multiple en- 
tangled pairs with fidelity above some threshold Fmin 
into a single entangled pair with higher purity||^||9|. The 
purification algorithm we use is described in detail in 
Refs.j^ |9|, 13^ |35|. For clarity we will present it in 
a form appropriate to the system under consideration, 
which uses repeated generation of electron spin entan- 
gled pairs to purify a stored nuclear spin entangled pair. 
Specifically, an electron spin entangled pair between sta- 
tions i and j is described by the density matrix diag- 
onal components {oe, 6e, Ce, de} in the Bell state basis 
{|vl/_),|$+),|<i>_),|vl/+)}, where 

l*±)!f = -^(|0a,)±|1.0,)) (30) 
\^±)'if = ^(|0.0,)±|ia,)). (31) 

We will refer to these diagonal elements as the "vector 
fidelity" J-e = {ce, 6e, Ce, de}, noting that the first el- 
ement {J^e)i = o,e encodes the fidelity with respect to 
the desired singlet state. A nuclear spin entangled pair 
between those stations is described by a similar vector 
fidelity J^n = {o-mbm Cn, (!„} in the nuclear Bell basis 



B. Purification 

To extend entanglement to long distances in the pres- 
ence of errors, active purification is required at each level 



l$±) 



in) 



1 

1 

71 



(li.T,) ± IT.i,» 
(IU,)±IT.T,». 



(32) 
(33) 
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The purification protocol calls for a local rotation of each 
spin system at both locations: 





- T2 - 

V 


^*|i)..) 


(34) 








(35) 


1 0-i^j 


- 71 (it).. - 


^*lt).,) 


(36) 


It)., - 


>7l(lt).^ 


-iU),j) , 


(37) 



followed by a two-qubit gate at each location: 



lii). 
ITO), 
IT 1). 



110) , 

111) . 

- IT 1), 

- IT 0), 



U 1), 
IT 0),. 
IT 1), 



IT i>, 

IT 0), 
ITi>,. 



(38) 



After these operations, the electron spin is projectively 
measured at both locations. When the two electron spins 
are in the opposite state, the purification step succeeds, 
mapping the remaining nuclear spins to a diagonal state 
Wn,b'„,c'^,d'„} with a; > a„. 

This purification protocol can correct any type of error, 
but it functions best for phase errors, which correspond 
to diagonal elements of the form {/, 0, 0, 1 — /}. To quan- 
tify the type of error associated with our entanglement 
generation scheme, we define a shape parameter v such 
that the vector fidelity for entangled spin pairs between 
adjacent nodes is 

To = {Fo, (1 - Fo)v, (1 - ^^0)^-, (1 - Fo){l - 2v)}. (39) 

Note that f — > corresponds to phase errors while 
y — > 1/3 corresponds to a Werner state with equal distri- 
bution of all error types. Note also that the assumption of 
diagonality imposes no restriction on the entangled states 
we generate, as any off-diagonal elements in their density 
matrices can be eliminated by performing random rota- 
tions (similar to the procedure for creating Werner states 
but without the synimetrization step) l35| . Furthermore, 
even without a randomization step, the average fidelity 
is determined by the diagonal elements |^. 



C. Errors 

In the presence of local errors in measurements and 
operations, the purification and swap procedures deviate 
from their ideal effect. To describe this we use the error 
model described in fsij . Measurement errors are quanti- 
fied using a parameter 77 such that measurement projects 
the system into the desired state with probability 77 and 
into the incorrect state with probability I — rj. For ex- 
ample, a projective measurement of state |0) would be 



Po^v\0) (0| + (1-7?)|1)(1| 



(40) 



Errors in local operations are accounted for in a similar 
manner. With some probability p, the correct operation 
is performed; otherwise one traces over the relevant de- 
grees of freedom in the density matrix and multiplies by 
the identity matrix (for further details seelsSj and refer- 
ences therein). For example, the action of a two qubit 
operation Uij would become 

U.,pUl ^ p U.,pUl + ^Tn, (p) ® (41) 

In our calculations, we neglect errors in single qubit op- 
erations and focus on two-qubit errors, which are likely 
to yield the dominant contribution. 

These errors determine the level of purification which 
is possible given infinitely many purification steps. They 
also determine how much the fidelity degrades during the 
entanglement swap procedure. Below we describe a re- 
peater protocol which, compared to the original proposal 
[3, reduces the required number of qubits at each re- 
peater station at the expense of extra connection steps. 
Owing to these extra connection steps, our protocol is 
slightly more sensitive to local errors than the original 
scheme. 



D. Nesting Scheme 

Previous proposals for fault tolerant long distance 
quantum communication have required larger and larger 
numbers of qubits at each node as the communication 
distance is increased. Here we describe a nesting scheme 
which can be used to communicate over arbitrarily long 
distances while maintaining a constant requirement of 
only two qubits per node. 

The scheme for nested entanglement purification is il- 
lustrated in Figure 2| For clarity, we will label puri- 
fied pairs by "A", pairs to be purified by "B", and auxil- 
iary pairs used to perform purification by "C" . Briefiy, 
an entangled pair ("B") is stored in the nuclear spins 
while an auxiliary entangled pair ("C") is generated in 
the electron spins. The purification protocol described 
ill iMj is then performed by entangling the electron 
and nuclear spins via the hyperfine interaction, and sub- 
sequently measuring the electron spins. Comparison of 
the measurement outcomes reveals whether the purifica- 
tion step was successful, resulting in a new stored pair B 
with higher fidelity. After successfully repeating the pro- 
cedure for m purification steps, (a technique sometimes 
referred to as "entanglement pumping" ) , the stored pair 
becomes a purified ("A") pair, which can then be used to 
create B and C pairs over longer distances. We may thus 
generate and purify entanglement to arbitrary distances. 
This procedure is analogous to the scheme in Ref. Q , but 
avoids the increase in the number of qubits required for 
that proposal. 

Mathematically, the scheme can be explained most eas- 
ily using inductive arguments. Suppose that we have 
a means to create and purify entanglement over k = 
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FIG. 4: Nesting scheme for generation and purification of en- 
tangled nuclear and electron spin pairs. In each node, the nu- 
clear spin degree of freedom is represented by the upper (black) 
circle, while the electron degree of freedom is represented by 
the lower (red) circle. Entanglement between different nodes 
is represented by a line connecting them. Ovals represent en- 
tanglement swap steps, and rectangles represent entanglement 
purification steps. For n — 2 the B and C pairs may be di- 
rectly generated. For n > 3, the first step illustrates how the 
B pair is generated, while the remaining two steps illustrate 
how the C pair is generated while storing the B pair. The 
arbitrary distance algorithm works for n > 6. 



2,3,- •• ,n/2 repeater stations {{n + l)/2 if n is odd), 
and that we know the vector fidelity J^A{k) and the time 
7^(/c) required for each distance. We can then deter- 
mine the time required and the vector fidehty possible 
after purification over n repeater stations. 

We begin by creating two purified nuclear spin A pairs 
over half the distance and connecting them via a central 
electron spin pair of vector fidelity J-q. In the presence of 
local errors, this yields a nuclear spin B pair with vector 
fidelity 



J^B{n) = C { {J^A(^),^o,-?^A(y)},?7,P 



(42) 



Here, C gives the vector fidelity obtained upon connecting 
the entangled pairs in the presence of local errors [s^, and 



n/2 and n'/2 are understood to represent (n — l)/2 and 
(n + l)/2 when n is odd. The B pair is created in a time 



X.(y)+To 



(43) 



where Tq is the time required to generate nearest- 
neighbour entanglement and tc is the classical commu- 
nication time between adjacent stations. We neglect the 
time required for local operations and measurement since 
these times are short compared to Tq and ntc. Similarly, 
we can find the vector fidelity and time for the electron 
spin C pair 



Tc{n) = r^(y-l)+To + (n-2)t,. 



(44) 



After performing one purification step, we obtain a nu- 
clear spin pair Ai , with vector fidelity determined by the 
purification function V 



(45) 



On average, the time required to perform this single step 
is 



{TB{n)+Tc{n) + {n-l%) 



(46) 



where Ps is the probability that the purification step 
succeeds. 

After m successful purification steps, the vector fidelity 
of the nuclear spin Am pair is 



^A,r,{n) =V{J^A,„^An),^c{n'),'n,p), 
and the average time required for its creation is 
{TA„^_,{n)+Tc{n) + {n-l%) 



Ps{m) 

(rc(n) + (n-l)t,)^n 
/ 1 



(47) 



(48) 



m m 



n=l k= 



1 



Ps{m) 



-rB{n)\{ 



k=l 



\Ps{m) 



(49) 



where Ps{m) — Ps{^A,n-i{n),Tc{n)). If we stop puri- 
fying at some fixed number M of purification steps, then 
the desired vector fidelity and time over distance n are 
given by 



TA{n) 
TAin) 



'^AmM. 



(50) 
(51) 



To complete the inductive argument, we must show 
that the protocol works for small distances. There are 
many schemes one can use to generate and purify en- 
tanglement over shorter distances, and one possibility is 
illustrated in Figure 01 In fact, once the physical param- 
eters for an implementation are determined, it should be 
possible to optimise the few-node scheme to minimise the 
required time or maximise the resulting fidelity. 
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E. Fixed point analysis 

As the number of purification steps increases to — > oo, 
the fidehty of the resulting entangled pair saturates. This 
saturation value can be found using a fixed point analysis 
(as described injs^l) by solving for the vector fidelity 
which is unchanged by further purification steps 

Ta=V{Ta.Tc.t^,V). (52) 

where we have explicitly included the local errors in the 
purification function V . This yields a fixed point fidelity 
^Fp{J'Ci V^P) which is independent o{Ta- Since the vec- 
tor fidelity J- a has three independent parameters char- 
acterising the diagonal elements of the density matrix, 
one might miss the fixed point. However, as the number 
of purification steps increases our simulations do indeed 
approach the calculated fixed point. We therefore calcu- 
late the fixed point as a function of distance to find the 
upper bound on the fidelity which can be attained for 
given Tq, L/Lq,p, and 77. 

F. Asymptotic Fidelity 

As the distance increases L — > 00, the fixed point fi- 
delity can approach an asymptotic value Foo- We can 
understand the existence of and its value by examin- 
ing the protocol as a function of nesting level. In partic- 
ular, to generate entanglement over n repeater stations 
we operate at nesting level i ^ log2 n, where we obtain a 
purified pair 

^i'^ =^FP(4'^^.p): (53) 

where J-pp is the fixed point solution to Eq. (|52|) . and 
{J^Fp)i — Fpp is the fixed point fidelity. We will then use 
this purified pair J-'^'^ to build up an auxiliary C pair on 
the next nesting level Since the fidelity over distance 
rt — 1 is greater than that over distance n, i.e. (Jvi(rt — 
l))i ^ i^Ai'n))i, the auxiliary pair fidelity we obtain will 
be greater than or equal to the first component of 

^ C{{J^^2\^a\^o,:Fo,J'o},V,p), (54) 

where C is again the connection function. This auxil- 
iary pair will then determine = J^Fp{TQ~^^\ri,p). 
When J^^2+'^ Tf , we have reached the asymptotic fi- 
delity = {^a)\ (see Figure E)), which is given by the 
intersection of the purification curve Eq. (|53|l and the 
auxiliary pair creation curve Eq. (|54|) . 

As was the case for the fixed point analysis, we must ac- 
count for all diagonal components of the density matrix in 
the Bell state basis (not just the fidelity a). Consequently 
the asymptotic fidelity represents an upper bound to 
which the system may converge in the manner indicated 
by our simulations. Finally, we should stress that our 
calculations have not incorporated loss due to the long 
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FIG. 5: Approach to asymptotic fidelity. The solid curve 
shows the purified fidelity obtained from the auxiliary pair, 
while the dotted curve corresponds to the auxiliary pair ( con- 
structed from two smaller purified pairs) on the next nesting 
level. The system moves between the curves at each nest- 
ing step, and the upper intercept of the two curves gives the 
asymptotic fidelity. For this calculation Fq = p = ri = 0.99, 
and V = 0. 

but finite memory time in the nuclear spins. This loss 
increases with the total time required for repeater opera- 
tion, and sets the upper limit on the distance over which 
our scheme could operate. 



G. Results 

The discussion of final fidelity may be summarized as 
follows: the fidelity obtained at the end of this nested 
purification procedure, F{m, L, Fo,p,r]), depends on the 
number of purification steps m, the distance L between 
the outer nodes, the initial fidelity Fq between adjacent 
nodes, and the reliability of measurements 77 < 1 and lo- 
cal two-qubit operations p < 1 required for entanglement 
purification and connection |3^ . As the number of pu- 
rification steps increases to ^ 00, the fidelity at a given 
distance L approaches a fixed point F Fpp{L, Fq,p, rj) 
at which additional purification steps yield no further 
benefit js^. Finally, as L increases, the fidelity may ap- 
proaches an asymptotic value Fpp — > Foo{FQ,p,rj). Fig- 
ure EJi illustrates the efhciency of the purification pro- 
tocol: for initial fidelities Fq > 97%, three purification 
steps suffice to produce entanglement at large distances. 

FigureElD demonstrates that our scheme permits gener- 
ation of high-fidelity, long distance entangled pairs in the 
presence of percent-level errors in polynomial time. Be- 
cause solid-state devices allow fast operations and mea- 
surements, the overall time scale is set by the classical 
communication time between nodes. As an example, us- 
ing a photon loss rate of ^ 0.2 dB/km and inter-node 
separation Lq ~ 20 km (so that in the limit of good 
detectors the collection efficiency is 10~°-^ ^ l/e)j a fi- 
delity set by an emission probability Pem ~ 8%, local 
errors rj = p = 0.5%, and just one purification step at 
each nesting level, our scheme could potentially produce 
entangled pairs with fidelity F ^ 0.8 sufficient to violate 
Bell's inequalities over 1000 km in a few seconds. For 



12 



0.94 
0.92 

10" 



i (a) 












" " ^FP . , 


)(v) 


20 40 60 80 100 
Distance (Lq) 


120 


(b) (V) 


(i) 










1 2 5 10 20 
Distance (Lg) 


50 100 



FIG. 6: (a) Fidelity scaling with distance. Points show results 
using 3 purification steps at each nesting level; dashed lines 
show the fixed point Fpp at each distance; dotted lines indicate 
the asymptotic fidelity Foo- For (a) and (b), measurements 
and local two-qubit operations rj = p contain 0.5% errors. The 
initial fidelity Fo is (i) 100% (ii) 99% (m) 98% (iv) 97% (v) 
96% with phase errors only, (b) Average time scaling with 
distance for m—3, given in units of To = (to+tc) /P , the time 
required to generate entanglement between nearest neighbors, 
and Lq, the distance between nearest neighbors. Measurement 
and local operation times are neglected. Note that the axes are 
logarithmic, so time scales polynomially with distance. 



comparison, under the same set of assumptions direct 
entanglement schemes would require ~ 10^" years. 

Fig- Ell shows that our scheme will operate in the pres- 
ence of 1 — p < 1% errors in local operations and percent- 
level phase errors in initial entanglement fidelity. Other 
types of error are in principle possible, and we consider 
nonzero shape parameters v for the initial fidelity J-q in 
Eq. (|39|l . The asymptotic fidelity shown in Fig. [7|3 indi- 
cates that, although the protocol we use is most effective 
for purifying phase errors, it also tolerates arbitrary er- 
rors. 



H. Optimization 

Once the parameters of the system are established, the 
protocol can be optimised to minimise the time required 
to generate some minimum fidelity Fmin over a distance 
L. We can vary the number of repeater stations ^ i/^o 
and the number of purification steps m (which need not 
be constant). We can also tailor the entanglement gen- 
eration procedure by changing the emission probability 
Pern to find the optimum balance between initial infi- 




Error in local operations 1-p Error in local operations 1-p 



FIG. 7; (a) Long-distance asymptote dependence on initial 
fidelity Fo of (i) 100% (ii) 99% (iii) 98% (iv) 97% (v) 96% 
with phase errors only, (b) Long-distance asymptote depen- 
dence on error type. For the calculations shown, Fq = 0.99, 
and the shape parameter ranges from v — Otov = 0.3. In both 
(a) and (b) measurement errors are set equal to operational 
errors, rj — p. 



delity 1 — Fq ~ Pem and entanglement generation time 
To cx 1/Pem- Finally, one could use more advanced opti- 
mal control techniques to vary the details of the protocol 
itself. In particular, it should be possible to speed up the 
algorithm by working simultaneously on multiple nest- 
ing levels, beginning entanglement generation and con- 
nection on the next nesting level as soon as the interior 
nodes are free. Further speed-up may also be possible 
in the case when collection efficiency is very high by us- 
ing coincidence detection in combination with e.g. time- 
bin encoding Q]. As noted previously such coincidence 
detection could also be advantageous for interferometric 
stability ,26., ..27„ .28.] . 

Ultimately, the speed of this protocol is limited by 
three factors: classical communication time between 
nodes, probabilistic entanglement generation, and se- 
quential purification. Faster techniques will require more 
efficient entanglement generation or larger numbers of 
qubits at each node to allow simultaneous purification 
steps. 



I. Comparison to other quantum repeater schemes 

This scheme combines the advantages of two pioneer- 
ing proposals for quantum repeaters 0, Early 
work showed that entanglement purification and swap- 
ping could be combined to permit efficient, fault-tolerant 
quantum communication over distances longer than an 
attenuation length This scheme incorporated error 
correction at the cost of increased physical resources, re- 
quiring nodes containing a number of qubits scaling at 
minimum logarithmically with distance ^3^. Owing to 
the difficulty of implementing even few-qubit quantum 
computation, implementation of this scheme remains a 
challenging goal. Our scheme is closely related to the 
original proposal with one key difference: by spatially 
rearranging the required physical resources, we can effi- 
ciently simulate their protocol while maintaining a con- 
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stant requirement on qubits per node. This makes our 
scheme amenable to reahstic physical implementation. 

Another physical implementation for quantum re- 
peaters uses atomic ensembles as a long-lived memory 
for photons p"5i| . Entanglement is generated by interfering 
Raman scattered light from two ensembles. The entan- 
glement is probabilistically swapped using an EIT read- 
out technique. This scheme elegantly avoids effects of 
the dominant photon loss error by conditioning success 
on photon detection. Our scheme primarily differs from 
this proposal in two ways: first, access to two-qubit oper- 
ations between electron and nuclear spin permits deter- 
ministic entanglement swapping; second, the two-qubit 
nodes allow active correction of arbitrary errors. 

V. PHYSICAL SYSTEMS 

We conclude with three specific examples for potential 
implementation of the presented method. 



A. Implementation with NV centers 

The NV center level structure illustrated in Fig. [21 al- 
lows implementation of all steps in the repeater protocol. 
The cycling transition from |0) to |e) is used for electron 
spin initialization by measurement, entanglement gener- 
ation, and electron spin state measurement. A series of 
ESR and NMR pulses can be used to perform arbitrary 
gat es between the electron spin and an adjacent ^■^C spin 
[33j. Consequently, nuclear spin state initialization and 
measurement is achieved by initializing the electron spin, 
mapping the nuclear spin state onto the electron spin, 
and subsequently measuring the electron spin. Entangle- 
ment propagation and purification can be implemented 
in NV centers by driving ESR and NMR transitions and 
using optical detection of the electron spin states. Once 
electron spin entanglement is established between nodes 
Ri and it can be transferred to the nuclear spins, 

leaving the electron degree of freedom free to generate en- 
tanglement between station Ri and Ri+i- Provided that 
we can reinitialize the electron spin without affecting the 
nuclear entanglement, we can perform the same proba- 
bilistic entanglement procedure. Note that ESR multi- 
plexing is required to perform a 7r/2 pulse independent 
of the nuclear spin; this can be accomplished simply by 
applying two ESR pulses at the two transition frequen- 
cies. 

We now consider the feasibility of implementing our 
repeater protocol using NV centers in diamond. Owing 
to the overlap of electron wavefunctions in the ground 
and excited states, most of the NV center optical emis- 
sion goes into the phonon sidebands. Other color centers 
in diamond, for example the NE8 center [36l Is^ may 
suffer less from this drawback. To enhance the relative 
strength of the zero-phonon line, it will be necessary to 
couple the NV center to a cavity. For NV centers coupled 



to cavities with Purcell factors ^ 10 we find that 
the dominant source of error is electron spin decoher- 
ence during the classical communication period. Using 
an emission probability Pem ^ 5%, a collection efficiency 
e ~ 0.2, and a classical communication time of tc ^ 70/iS 
over Lq ^ 20 km, we find the fidelity of directly entan- 
gled pairs can reach Fq ^ 97% for electron spin coherence 
times in the range of a few milliseconds. Electron spin co- 
herence times in the range of 100/xs have been observed at 
room temperature and significant improvements are ex- 
pected for high purity samples at low temperatures |38j |. 
The large hyperfine splitting allows fast local operations 
between electron and nuclear spin degrees of freedom on 
a timescale ^ 100ns 33] much shorter than the decoher- 
ence time, allowing l—p< 1%. Finally, cavity enhanced 
collection should significantly improve observed measure- 
ment efficiencies of ry ~ 80% [3^ . 



B. Alternative implementation: quantum dots 

Our discussion thus far has attempted to remain gen- 
eral while exemplifying our proposal using NV centers. 
The basic idea of using two-qubit repeater stations should 
be applicable to a wide variety of systems featuring cou- 
pled electron and nuclear spins. To illustrate an alterna- 
tive implementation, wc consider doped self-assembled 
quantum dots whose electron spin is coupled to col- 
lective nuclear states in the lattice. Compared to NV 
centers, this system offers large oscillator strengths and 
the potential for Raman manipulation. Doped semicon- 
ductor quantum dots have been considered in a vari- 
ety of quantum computing proposals and related tech- 
nologies |3^ ^3 ■ The spin state of the dopant electron 
provides a natural qubit with relatively long coherence 
times. Assumin g a high degree of nuclear spin polariza- 
tion (Fn > 0.95) and active ESR pulse correction, the 
electron spin dephasing time is expected to be 1 ms '4^ . 
The spins of lattice nuclei in the quantum dot provide 
an additional, quasi-bosonic degree of freedom with ex- 
tremely long coherence times ('^ 1 s with active correc- 
tion [43j ) . Such ensembles of nuclear spin have been con- 
sidered for use as a quantum memory |44| and, by taking 
advantage of the non-linearity of the Jaynes-Cummings 
Hamiltonian, as a fundamental qubit for a quantum com- 
puter 45j. 

Unlike the spin triplet state of the NV centers, the 
conduction band electron has two states, |t) and ||), cor- 
responding to spin aligned and anti-aligned with an ex- 
ternal magnetic field Bext\\z. The quantum dot system 
also differs from NV centers in that it can be manipu- 
lated using Raman transitions: when the external field 
and growth direction are perpendicular (Voigt geometry) , 
two allowed optical transitions to a trion state produce a 
lambda system ; moving towards aligned field and growth 
directions (Faraday geometry) suppresses the "forbidden" 
transitions, as shown in Fig. 

Electron spin coherence can thus be prepared via Ra- 
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FIG. 8: (a) Level structure for single electron to trion transi- 
tion in a single- electron doped HI- V or II- VI quantum dot with 
external magnetic field in close to a Faraday geometry with 
Zeeman splitting ujz, heavy-hole splitting uj^ and up to four 
optical fields of different frequency and polarization. Dashed 
lines indicate weak dipole moments due to small magnetic field 
mixing. The triplet two-electron states are not included due 
to a large (> ImeV) exchange energy allowing for complete 
suppression of their effects, (h) Electronic (\'\) and col- 
lective nuclear states and their transitions for singly 
doped quantum dots in a polarized nuclear spin lattice. 



man transitions or by standard ESR setups, and changes 
in effective magnetic field can be accomplished by off- 
resonant, spin-dependent AC Stark shifts with ctj^ light. 

Although optical transitions in doped quantum dots 
can exhibit homogeneous broadening T ^ 100 GHz~ 
10— IOO7 ^4^, the corresponding error can be made neg- 
ligible by sending the output from the cavity through a 
frequency filter with a linewidth of a few hundred MHz. 
[sojl Moreover, we note that InAs quantum dots have been 
successfully coupled to microcavities with Purcell factors 

~ 10 im. 

Whereas the NV center electron spin was coupled to 
a single nuclear impurity, the electron in a quantum dot 
couples to collective excitations of many thousands of 
nuclei. We briefly discuss this system; further details are 
given m Refs. lillill. The Hamiltonian governing this 
interaction is 



k k 

where 7^ is the gyromagnetic ratio for nuclear spin / , the 
nuclear spin coupling amplitudes satisfy oc | Pi 
Sfc -^fc ~ 1' ^"^^ ^ ~ A/fiJ^k^k (A- is the hyperfine 
interaction constant). By identifying collective nuclear 

spin operators, A = J^k ^kl, the hyperfine term may be 

written MIS ■ A. For simplicity we restrict the following 
discussion to the case of perfect nuclear polarization (so 
the initial state of all Nn nuclear spins in the quantum 
dot is |0) = |-/) (g) ... (g) |-/) for /-spins). Then an 
effective Jaynes-Cummings type Hamiltonian describes 
the system: 



Hl^f = hujl^f + hn/2[A+S- + A^S+). 



(56) 



with corrections of order fi/ Afc ~ A/Nn. The effec- 



tive Zeeman splitting oj 



elf 



■ lA/h is dominated by 



the field associated with the polarized nuclear spins ; for 
example, in GaAs quantum dots, this Overhauser shift 
is lA/h ^ 33GHz. The large detuning u^^f suppresses 
interactions which exchange energy between the electron 
and nuclear spins. By changing the effective magnetic 
field, we can shift the system into resonance uj^-^^ — > 
to drive Rabi oscillations between the electron spin and 
the collective nuclear state, see figure IHb. Pulsing the 
appropriate effective field permits a controllable map be- 
tween the electron spin state and the collective nuclear 
degrees of freedom spanned by |0) and |1) = A+ |0). In 
addition, more complicated sequences of electron-nuclear 
spin interaction and electron spin manipulation allow for 
arbit rary two qubit operations on |t) , |i) and |0) , |1) (see 
Ref. |45j | and commentary therein). 

Measurement and initialization proceed in much the 
same manner as described for NV centers. The state of 
the electron spin system can be read out by exciting a 
cycling transition with resonant cr+ light, and measure- 
ment and ESR (or Raman transitions) can be employed 
to initialize the system in the desired state. As the effec- 
tive Knight shift of the electron spin is negligible on the 
time scales of entanglement preparation, the collective 
nuclear state's coherence is unaffected by this process. 
Due to the improved selection rules and possibility of 
Raman transitions, it may be more effective to use the 
Raman entanglement generation scheme. 

The nuclear state of the quantum dot can be prepared 
by cooling the nuclear spins using preparation of electron 
spin and manipulation of the effective magnetic field |48j | . 
In practice, this leaves the nuclear system in a state I'D) 
with the same symmetry properties as the state |0) de- 
scribed above j42,J. To date, 60% nuclear spin polar- 
ization has been achieved by optical pumping in GaAs 
quantum dots As was the case with NV centers, the 
nuclear spin state can be read out by preparing the elec- 
tron spin in the ||) state, mapping the nuclear state to 
the electron spin state, and measuring the electron spin 
state. 



C. Atomic Physics Implementation 

Compared to the solid state implementations we have 
considered so far, implementations in single trapped 
atoms or ions have the advantage that they typically have 
very little broadening of the optical transitions. Because 
atomic systems do not reside in a complicated many-body 
environment, their internal degrees of freedom can have 
very long coherence times. For most atomic systems, 
however, it is hard to identify a mechanism which allows 
one degree of freedom, e.g., the nuclear spin, to be decou- 
pled while we probe some other degree of freedom, e.g., 
the electron spin. Below we describe a system which does 
fulfill this requirement, although practical considerations 
indicate implementation may be challenging. 

We consider alkali-earth atoms, such as neutral mag- 
nesium, and chose an isotope with non- vanishing nuclear 
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spin (^Mg). The lowest lying states of magnesium are 
shown in Fig. |51 (electronic structure only) . Instead of 
the electronic spin states we have considered so far, i.e., 
for NV centers and quantum dots, we will use states 
which differ both in spin and orbital angular momen- 
tum. The stable ground state ^So will serve as state |0). 
In this state, the electronic degrees of freedom have nei- 
ther spin nor orbital angular momentum and the nuclear 
spin is thus decoupled from the electronic state. The ex- 
cited state "^Pq (whose hyperfine interactions also vanish 
to leading order) will provide state |1). Note that the 
triplet-singlet transition from ^Pq to the ground state 
is highly forbidden and this state has an extremely long 
lifetime, but transitions between the two states can still 
be induced with a strong laser. 

To create entanglement we couple the ground state to 
the excited state ^Pi with a laser field and collect the 
scattered light. From this excited state the atom es- 
sentially always decays back into the ground state. If 
the driving is detuned much further than the hyperfine 
splitting in the excited state, the nuclear spin is also de- 
coupled during this process. The nuclear spin can there- 
fore by used to store information while we entangle the 
electronic state with another atom. Finally, to imple- 
ment gates between the electronic and nuclear states one 
should, for instance, couple the |0) state to another state 
in the atom where there is a hyperfine interaction, for 
example using resonant excitation of the ^Pi state. 



s=o 



(3s3p)^P, 



(3S^)'S„ 



|0) 



S=l 



-J=2 

-J=l (3s3p)'p; 
-J=0 



FIG. 9: Electronic level structure of atomic magnesium. The 
electronic ground state ^So has vanishing spin and orbital an- 
gular momentum. In this state the nuclear spin therefore de- 
couples from the electronic degrees of freedom. An entangling 
operation which is also insensitive to the nuclear degree of 
freedom can be achieved with lasers which are detuned much 
further than the hyperfine splitting in the ^Pi level. 

Finally, we note that all three physical implementa- 
tions we suggest operate in the visible or near-IR, and 
will likely require high-efficiency frequency conversion to 
telecom wavelengths for low-loss photon transmission. 



VI. CONCLUSION 

In conclusion, we propose a method for fault tolerant 
quantum communication over long distances requiring 
only probabilistic nearest-neighbor entanglement gener- 



ation, two-qubits per node, and two-qubit operations. 
We compare several schemes for entanglement genera- 
tion and discuss two solid-state systems and an atomic 
system which might be used to implement them. Po- 
tential applications include secure transmission of secret 
messages over intercontinental distances. 
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VII. APPENDIX 

Photo-bleaching is a detectable error, so it does not 
affect the fidelity of entanglement generation or mea- 
surement, as we described above. However, it can in- 
crease the time required for these operations. Fluores- 
cence correlation experiments are consistent with assign- 
ing a metastable singlet structure to the shelving state, 
which is coupled strongly to the Mg = ±1 excited states 
but only weakly to the Mg — excited state (see Fig- 
ure [21), |2^]. We need to account for the possibility that 
our NV center bleaches during entanglement generation, 
requiring us to start over. During each attempt we reso- 
nantly excite the Mg = transition with some probabil- 
ity Pcm (see Eq. El . To quantify the population lost to 
\W) we will consider a model where the Mg = ±1 excited 
states decay to the shelving state at rate Fg. The oscil- 
lator strength for the Mg = ±1 optical transitions are 
unknown, so we will assume that the Rabi frequencies on 
the Mg — ±1 transitions are fl' . 

During one attempt at entanglement generation, the 
probability to end up in the shelving state is 



w 



(57) 



where 5 is the detuning from the Mg = ±1 optical tran- 
sition. (The excited state energies are strongly inhomo- 
geneously broadened, so 5 is not precisely known; this 
detuning should be controllable using strain or applied 
electric fields.) On average, a large number of attempts 
are required for successful entanglement gen- 
eration. Consequently the total probability for the sys- 
tem to end up in the shelving state during entanglement 
generation is 



Pi 



w 



AVs- 



7 



52 ^2^2 



(58) 



where /i(/i') is the oscillator strength for the Mg — 0(±1) 
transition. The precise values of these parameters are 
unknown, but we can estimate their order of magnitude: 
Fs - 1 - 10 MHz, 7 -t- F ~ 100 MHz, 6 ~ 1 GHz, /z' - /x, 
yielding 



Pw 



4(10- 



10- 



(59) 
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If this error rate is too large, we can also check for photo- 
bleaching at intervals during the entanglement proce- 
dure. 

The shelving state poses a similar problem during mea- 
surement. In this case, the Mg = transition is strongly 
illuminated so that at least one photon reaches the de- 
tectors: Pem ~ Under the same illumination, any 
population in |1) will end up in the shelving state with 
probability 

^ ^ v^n'^ r7 m'^ (10-3-10-4) , ^ 



Note that the measurement fidelity is unaffected by 
photo-bleaching if we verify that the center is optically 
active by observing fluorescence either directly from the 
|0) state or after applying a multiplexed ESR pulse to 
the |1) states. Ultimately, in this model the only effect of 
the shelving state is to reduce the success rate for entan- 
glement generation and measurement by of order a few 
percent. Finally, we should note that the eff'ect of the 
shelving state on the nuclear spin state is currently not 
known, and could potentially complicate the sequence of 
operations necessary upon detection of a shelving event. 
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